Non-exponential London penetration depth in RFeAsOo.gFo.i (R=La,Nd) single 
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The superconducting penetration depth, X(T), has been measured in RFeAsOo.gFo.i (R=La,Nd) 
single crystals (R-llll). In Nd-1111, we find an upturn in A(T) upon cooling and attribute it to 
the paramagnetism of the Nd ions, similar to the case of the electron-doped cuprate Nd-Ce-Cu-O. 
After the correction for paramagnetism, the London penetration depth variation is found to follow a 
power-law behavior, AA_l(T) oc T 2 at low temperatures. The same T 2 variation of A(T) was found 
in non-magnetic La-1111 crystals. Analysis of the superfluid density and of penetration depth 
anisotropy over the full temperature range is consistent with two-gap superconductivity. Based on 
this and on our previous work, we conclude that both the RFeAsO (1111) and BaFe2As2 (122) 
families of pnictide superconductors exhibit unconventional two-gap superconductivity. 
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PACS numbers: 74.25.Nf,74.20.Rp,74.20.Mn 

A year after the discovery of iron pnictide superconduc- 
tors [l| , the order parameter (OP) symmetry still remains 
one of the most important open questions. The pnictides 
are complex superconductors with a possible influence 
from magnetism, similar to the cuprates |2j, and possi- 
ble multigap superconductivity, as observed in MgB2 Q • 
Theoretical reports predict that a spin-density wave in- 
stability and multiple Fermi surface sheets can give rise to 
a complex momentum dependence of the superconduct- 
ing OP. It can have opposite signs on different sheets of 
the Fermi surface (FS) and have s-wave (s ± -wave) J], 
anisotropic s-wave [5], extended s-wave with nodes 6|, 
or dj.2_ a 2-wave symmetry Q. 

The London penetration depth, Al(T), is among the 
most useful probes of the OP symmetry A fully 
gapped FS leads to exponential saturation of \l(T — > 0). 
The saturation temperature is determined by the mini- 
mum value of the gap, which can either be due to pairing 
anisotropy or/and different gap amplitudes on different 
sheets of the FS. The presence of the gap zeroes (point 
or line nodes) on the FS leads to a non-exponential vari- 
ation of Al(T), e.g. power-law behavior. In the high- 
T c cuprates, precision measurements of \l(T) were the 
first to convincingly show the existence of rf-wave pair- 
ing state Measurements of A(T) in recently dis- 
covered FeAs-compounds reveal a contradictory picture. 
Extensive studies by our group on both electron-doped 
Ba(Fe,Co) 2 As 2 10] and hole-doped (Ba,K)Fe 2 As 2 
(122) systems have found Al(T) oc T 2 behavior down to 
0.02T C , which may be a signature of unconventional 
12, |l3|. Microwave measurements in hole- 
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doped 122 suggest a fully gapped state with two gaps 
|14j . whereas /xSR measurements show a linear temper- 
ature dependence of the superfluid density, suggesting a 
nodal gap [Hj|. In R-llll, an exponential dependence 
in A(T) was reported in single crystals of Sm-1111 [Til ] 
and Pr-1111 13, LL8| as well as in our own experiments 
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FIG. 1: (Color online) Transverse scans through the (200) 
and (220) Bragg reflections of the Nd-1111 single crystal. The 
lines are hts to Lorentzians yielding a full width half maxi- 
mum of 0.019 deg and 0.028 deg, respectively. Left inset: 
photographs of Nd-1111 (top) and La-1111 (bottom) single 
crystals. Right inset: A(T) in the full temperature range for 
three samples. 



on Nd-1111 Il9l . A nodeless gap in R-llll is supported 
by ARPES [20( and by point contact spectroscopy with 
sing le 21 1 and two-gaps @, whereas [iSR 23] and NQR 
24( results could be interpreted as both, multigap s- or 
d-wave state. Furthermore, penetration depth measure- 
ments in stoichiometric 1111 LaFePO suggest the exis- 
tence of line nodes [25 1. 

In this paper we report on penetration depth measure- 
ments in RFeAsOo.gFo.i (R=La, Nd) single crystals. The 
Nd-1111 crystals were much larger than those used in our 
previous study (l9j and therefore our new data were ob- 
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tained with notably improved signal-to- noise (S/N) ra- 
tio. We find that the measured A(T) is not exponen- 
tially flat, but shows a paramagnetic upturn caused by 
the Nd 3+ ions. Similar situation was resolved ten years 
ago in the electron doped cuprate Ndi.ssCeo.isCuO^ 

(ncco) r 



271 ] . Correcting the data for this para- 



magnetic contribution, we have determined the Lon- 
don penetration depth, \l(T), which does not saturate 
down to 0.02T C and is best described by the power-law 
AAl(T) oc T 2 . Measurements on the recently grown 
non-magnetic La-1111 crystals reveal the same power- 
law without any corrections. The anisotropy of the pen- 
etration depth, 7a = A c /A a b, was found to increase upon 



cooling, similar to the 122 compounds [11|, |28[, which 
combined with the positive curvature of the superfluid 
density at elevated temperatures may be considered a 
consequence of multiband superconductivity. 

Single crystals of RFeAsOo.gFo.i (nominal O and F 
composition) were extracted from 5 mm diameter pellets 
obtained using high pressure synthesis [2^]. Individual 
superconducting single crystals have dimensions up to 
650 x 180 x 120 fim 3 for Nd-1111 and 330 x 240 x 10 fim 3 
for La-1111 (left inset of Fig. [T|). In the following, we 
present results for the largest samples with the highest 
S/N ratio. Similar results were obtained in four other 
somewhat smaller crystals. The mean T c was 45 K for 
Nd-1111 and 14 K for La-1111 crystals. 

The samples were characterized using synchrotron X- 
rays (6ID-D beamline in the MUCAT sector at the 
Advanced Photon Source, Argonne) with an energy of 
99.6 keV and an absorption length of 1.5 mm, probing 
throughout the sample thickness. The incident beam was 
collimatcd to 0.1x0.1 mm 2 and the sample c-direction 
was aligned parallel to the beam. Two-dimensional scat- 
tering patterns were measured with a MAR345 image- 
plate positioned 1705 mm behind the sample. Complete 
reciprocal planes were recorded by rocking the sample 
during the exposure through two angles perpendicular to 
the beam over a range of ± 2.4 deg [3C|. This way the 
Bragg reflections within reciprocal planes perpendicular 
to the incoming beam were recorded in a single expo- 
sure of the detector (typically 3.5 min). Stepwise trans- 
lations of the sample perpendicular to the incident beam 
allowed measurements of the spatially resolved diffrac- 
tion patterns in the basal ab-plane. The entire sam- 
ple was scanned and no impurities, misoriented grains or 
diffuse signals from disordered material were detected. 
The excellent quality of the single crystals used in this 
study can be seen from transverse scans, at room temper- 
ature, through selected Bragg reflections extracted from 
the two-dimensional pattern, Fig. Q] We find sharp peaks 
with widths of 0.019 deg and 0.028 deg, typical of high 
quality single crystals. The diffraction pattern did not 
change qualitatively on cooling down to 15 K, revealing 
no traces of the structural transition that is observed in 
the non-superconducting parent compound. 




FIG. 2: (Color online) (a) Frequency shift A/(T) for a Nd- 
1111 crystal (symbols), contribution from the sample holder 
(continuous line) and ^ffi (dashed line). Inset: A/(T) for 
the La-1111 crystal, (b) Low-temperature region of A(T) in 
the Nd-1111 crystal. The lines are fits to Eq. (fTa)) (dashed) 
and Eq. (|lb[) (continuous) Inset: Original data (symbols) and 
Al(T) obtained after dividing by ^/JL from a fit to Eq. (|la[) 
(dots) and to Eq. (|lb|) (continuous line). 



The magnetic penetration depth, A(T), was measured 
by placing the sample inside a 1 fiH inductor of a self- 
resonating tunnel-diode resonator (TDR) with resonant 
frequency f = l/2ny/LC w 14 MHz. The excita- 
tion ac magnetic field, H ac ~ 10 mOe, is much smaller 
than the lower critical field H c \ ~ 100 Oe, assuring 
that the sample is in the Meissner state. The mea- 
sured quantity is the shift in the resonant frequency 
A/ = j(T) - f = -4vrx(T)G, where X is the total 
magnetic susceptibility and G ~ foV s /2V c (1 — N) is a 
geometric calibration factor defined by the coil, V c , and 
the sample, V s , volumes and the demagnetization factor 
N. G is measured directly by pulling the sample out 
of the coil at the lowest temperature [8J . The suscepti- 
bility in the Meissner state can be written in terms of 
Aj,(T), /i(T), and a characteristic sample dimension w 

as 4tt X (T) = [Vm(T)A l (T)/ W ] tanh [^JIw/X L {T)] - 1, 

where fi{T) is the normal state paramagnetic permeabil- 
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ity |a lay- 

Figure [Ha) shows the frequency shifts for Nd-1111 
(main frame) and La-1111 (inset) crystals at low tem- 
peratures (the whole range is shown in the inset in 
FigJTJ). The samples were mounted with H ac \\ c, so that 
A/ oc AX a b- Note that the total frequency change over 
15 K is about 1 Hz, which is less than 0.1 ppm. However, 
this is still significantly larger than the noise level of our 
system, shown in Fig. [21(a) for comparison. Whereas the 
data in La-1111 are monotonic, the measured Af(T) in 
Nd-1111 shows an upturn below 4K. A similar upturn 
in the electron-doped cuprate NCCO was explained by 
the local moment magnetism of the Nd ions 26], 2t|. In 
the presence of paramagnetic ions, the screening length 
is reduced to X(T,fi) — Al(T)/ \J f^T), where Xl is the 
London penetration depth. At low temperatures, where 
tanh(w/A) rs 1, the measured Af can be written as 
Af{T) = GyfJlX L (T) [HI. If we assume Currie- Weiss be- 
havior for the magnetic susceptibility %(T) = C/(T+9), 
with C being the Curie constant (An factor included) 
and 9 being the Curie- Weiss temperature, we obtain an 
upturn in Af(T) provided the paramagnetic increase is 
larger than the decrease due to Al(T). Figure [21(a) shows 
that ^ fi(T) clearly dominates the experimental data be- 
low 4 K. To subtract the contribution of the paramag- 
netic term, we fit the data assuming Xl(T) to be either 
exponential, Eq. ([la]) , or power law, Eq. ([lb"]). 
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A(T) = vV(T)A(0) [1 + AT n ] (lb) 

As shown in Fig. [U(b) , f° r both forms the best fits are 
obtained for 9 0.2 K, which suggests the possibility of 
antiferromagnetic order of the Nd moments. A fit to the 
s-wave BCS form Eq. (fTa|) with Ao and the Curie- Weiss 
constant C as the free parameters yields Ao = (0.8±0.2)T c 
and C=0.085±0.03. Through many iterations we find 
that the best convergence to Eq. (|lb[) is obtained for 
n 2. This power not only reproduces the data ac- 
curately, but it gives a narrow range of the variation for 
C = 0.18±0.02, which is close to the theoretical C = 0.25 
estimated for the free Nd 3+ moment of 3.62/is. The inset 
of Fig. [H(b) shows the London penetration depth Xl (T) 
obtained after dividing by \J /i(T) for two sets of values 
(C, 9): (0.085, 0.2) from the best fit to Eq. dTa]) and 
(0.18, 0.2) from the best fit to Eq. (Ib|). We notice that 
an assumption of exponential Al(T) does not remove the 
low temperature upturn. A(T) obtained from the fit to 
Eq. (fib)) docs not saturate exponentially for T < 6K, 
contrary to earlier reports that analyzed the total mea- 
sured penetration depth [ill [13, El • As shown in Fig. [3l 
Xl{T) is well described by a quadratic temperature de- 
pendence down to 0.02 T c for all studied R-llll crystals. 
This remarkable finding together with the results from 




0.00 



0.02 



0.04 



0.06 



0.08 



0.10 



(T/T c 



FIG. 3: (Color online) AX L vs (T/T c ) 2 after dividing the 
paramagnetic contribution for two Nd-1111 samples and raw 
data for a La-1111 crystal (continuous lines). The curves have 
been shifted vertically for clarity. 




FIG. 4: (Color online) The superfluid density, p(T), for R- 
1111 crystals calculated with the two values of 200 and 350 
nm for Al(0) = (symbols). Solid lines are fits to the two-gap 
model described in the text. Inset: temperature dependent 
anisotropy of the penetration depth. 



Ref. [HI EH, suggest that the nearly quadratic temper- 
ature dependence of the London penetration depth is a 
universal characteristic of the iron-pnictide superconduc- 
tors, not only of the FeAs-122 family. This finding offers 
compelling evidence for unconventional superconductiv- 
ity in RFeAsO . 9 F .i. 

Whereas the pairing mechanism in FeAs superconduc- 
tors remains an open question, we show below that the 
penetration depth anisotropy and superfluid density are 
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consistent with the existence of two gaps in R-llll. To 
calculate the anisotropy j\(T) — \ c /\ a b (both quantities 
are now London penetration depths after accounting for 
the paramagnetic contribution) , we need to know the ab- 
solute values of \ a b{T) and A C (T). For the in-plane pen- 
etration depth, /iSR experiments give A a &(0) = 200 nm 
|32l ]. To obtain the absolute value for A C (T), we use the 
anisotropy of H c2 (T) near T c as described in Ref. [33j |. 
From our TDR measurements 11911 as well as from from 
resistivity [34j | and specific heat [35j measurements, 7a = 
ld T < T c) « 4-5 at T c . The inset of Fig. g] shows ~f\(T) 
calculated assuming j\(T) = 4.5. Regardless of the ini- 
tial value chosen, jxiT) increases upon cooling, reaching 
7a (0) « 17 ±3 at T = 0. Such behavior of the anisotropy 
was first reported in Sm- and Nd-1111 and interpreted 
in terms of multiband superconductivity in Ref. [361 ] . A 
similar temperature dependence was also observed in the 
FeAs-122 superconductors j28|. A possible explanation 
for the increase of 7a (T) upon cooling is the existence of 
two superconducting gaps with different mag nitudes on 
the FS with different anisotropies 11 , 36 



While the low-temperature behavior of Ax, (T) is clearly 
non-exponential, at intermediate temperatures the re- 
sponse is dominated by the quasiparticles thermally ex- 
cited over the gaps. One can try a simple s-wave model 
to probe possible two-gap superconductivity, prompted 
by the positive curvature of p(T), similar to MgB2 
39]. Figure 0] shows the in-plane superfiuid density, 
p = A 2 (0)/A 2 (T). Due to uncertainty in A(0) we have 
used the two values of 200 nm and 350 nm that agree 
with the literature for the 1111 pnictides [23|, [32|. To 
fit the data, we have used the model, where p(T) — 
epi(Ai) + (1 — e) P2 (A2) with two clean s-wave gaps 
Ai (T) and A2 (T) along with relative densities of states, e 
and (1 — e). Successful fits were obtained in all cases with 
the following values of (Ai/fc^, A 2 //cb and e): La-1111 
(11.29 K, 4.16 K,0.83) for A(0) = 200 nm and (14.98 K, 
5.24 K, 0.82) for A(0) = 350 nm. For Nd-1111, we obtain 
(47.25 K, 14.71 K, 0.86) and (64.15 K, 21.61 K, 0.85). 
In all samples, the ratio between the gaps is A1/A2 « 3, 
with the large gap at the FS sheet with 80% of the total 
density of states. This result is similar to previous reports 
on Sm-1111 and Pr-1111 0, H Si] . This analysis 
poses two interesting questions. First, for all cases, we 
obtain Ai/fcsT c and A2 /ksT c smaller than the weak- 
coupling value of 1.76. Second, for the experimentally 
determined value of A a &(0)=200 nm, the superfiuid den- 
sity shows a pronounced positive curvature at interme- 
diate temperatures and a strong suppression near T c . In 
the vicinity of T c , the superfiuid density should be linear 
with temperature, p(T) = rj (1 — T/T c ). For a clean s- 
wave superconductor -q — 2, whereas for RFeAsOo.gF .i 
we obtain 77 ps 1. This reduced slope agrees qualita- 
tively with the theoretical calculations for the case of 
two s-wave gaps (40| and/or extended s + symmetry 13 1 



rj oc I/T12, therefore a higher t\i will lower the slope of 
p(T) at T c . 

In conclusion, we find that the in-plane London 
penetration depth in single crystals of RFeAsOo.gFo.i 
(R=La,Nd) follows a power-law temperature depen- 
dence, AA q6 (T) cx T 2 for T < T c /3. The penetration 
depth anisotropy 7a increases upon cooling and the anal- 
ysis of the superfiuid density suggests two-gap super- 



conductivity. Given our results on 122 family [UJ, lllj . 
we conclude that iron-based pnictides are unconventional 
two-gap superconductors. 
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when a large interband scattering rate T12 is considered: 
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